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^ | Chiral restoration at finite temperatures is studied in chiral quark models with 

non-local regulators. At the leading-iV c level we find transition temperature of the 
. order lOOMeV. Meson-loop contributions are also analyzed and found to have a very 

! small effect. 

in 

Recently, considerable activity has been focused on chiral quark models with 
non-local interactions Qg^ggg^ggJigjT^Q^Q^^^ . In this approach 



OS 

interactions depend on the momenta carried by the quarks, which leads to a 
momentum-dependent quark mass, generated by spontaneous breaking of the 
chiral symmetry. There are several important reasons why it is worthwhile to 
consider such models: (1) Non-locality arises naturally in the most sound and 
successful approaches to low-energy quark dynamics, namely the instanton- 
liquid model and the Schwinger- Dyson resummation techniques |2T 



^ ■ Hence, we should deal with non-local regulators from the outset and models 

with local interactions (the Nambu-Jona-Lasinio model and its progeny) may 
be viewed as simplifications to the non-local models. (2) The non-local regu- 
larization leads to preservation of anomalies |||22| and to correct charge quan- 



tization Hl5yi6[ . With local methods, such as the proper-time regularization 
or the quark-loop momentum cut-off P,|ID],^31B^] the matching of anomalies 
can only be achieved if the finite (anomalous) part of the Euclidean action is 
left unregularized, and the infinite (non-anomalous) part is regularized. When 
both parts are regularized, anomalies are violated substantially |25|.|2(J and 



the theory is not consistent. The different treatment of the anomalous and 
non-anomalous parts of the action is rather artificial and it is quite appealing 
that with non-local regulators both parts of the action can be treated in a 
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unified way. (3) With non-local interactions the effective action is finite to 
all orders in the loop expansion (i.e. the 1/N C expansion). In particular, the 
meson loops are finite and it is not necessary to introduce additional cut-offs, 
as was the case of local models |27],^,[2{|]. As a result, the non-local models 
have more predictive power. (4) In the baryonic sector the non-local models 
are capable of producing stable solitons |15| without the extra constraint that 
forces the a and tt fields to lie on the chiral circle. Such a constraint is external 
to the known derivations of effective chiral quark models. 

In view of the above-listed advantages of the non-local models over the con- 
ventional Nambu-Jona-Lasinio-like models it is desired to check other appli- 
cations of non-local effective chiral theories. In this paper we investigate the 
restoration of chiral symmetry at finite temperatures and vanishing baryon 
density. The basic quantity of our study is the quark condensate, which in the 
chiral limit of the vanishing current quark mass is the order parameter of chiral 
symmetry. We show that chiral restoration is a second-order phase transition, 
as in the local models, and it occurs at temperatures T c ~ lOOMeV, which is 
somewhat less than in local theories, where the typical T c is around 150MeV. 
We also investigate the dependence of the critical temperature T c on the de- 
tailed shape of the non-local regulator and find it to be weak. With a finite 
current quark mass the chiral transition is, as usual, a smooth cross-over. In 
addition, we show that the effects of meson-loops (higher-order effects in the 
1/N C expansion) of the quark condensate are very small in non-local models. 

We begin with a brief review of the model. Our starting point is the two-flavor 
effective non-local Euclidean action in the bosonized form 

/ = -Trln(-i7^ + m + r$r 4 ,r) + 4/rf 4 x($g + <l> 2 ), (1) 

where -j has the meaning of the coupling constant of the four-quark interac- 
tions, m is the current quark mass, $ = ($cb ®a) are the mean meson fields 
with quantum numbers of the a meson and the pions, respectively. The fields 
$ are local in coordinate space. Matrices T$ = (r$ ,r$ a ) are acting in spin 
and flavor spaces, with r$ = 1 and T$ a = i^T a . The regulator operator r is 
local in momentum space, i.e. (p\r\p') = 5(p — p')r(p). It introduces a cut-off 
scale A, limiting appropriately momenta of the quarks. The symbol Tr denotes 
the full trace, i.e. the functional trace and the matrix traces over the Dirac, 
flavor and color degrees of freedom. 

Through the use of the Euler-Lagrange equations for the $ field it is straight- 
forward to show that the action (0) is equivalent to the model where the 
four-quark interaction has a separable form namely 

J_ r d 4 px r d 4 p 2 r d 4 p 3 r d 4 p A 
int " 2a 2 J (2tt) 4 J (2vr) 4 J (2vr) 4 J (2tt) 4 X ^ 
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<5 4 0i +P2 +P3 +P4)^(Pi)^(Pi)r$r(p 2 )^(P2)^(P3)^(P3)r$r(p 4 )V'(p4). 



This form has a very clear interpretation via Feynman diagrams: the quark 
vertices carry a regulator on each quark line entering the vertex. 

In the instanton-liquid model the regulator can be evaluated for low Euclidean 
momenta. It has the form fllS. ISpU 



r(p 2 ) =-~(I (z)K (z)-h(z)K 1 (z)), z = ^, (3) 

where I and K are the modified Bessel functions, and p is the average instanton 
size, of the order of (600MeV) _1 . We introduce A = 2/p, which is therefore 
of the order of 1.2GeV. We will also study an approximate formula for the 
regulator given by the Gaussian function: 

r (f) = exp (4) 

We first analyze the model at the one-quark-loop level. In the vacuum the val- 
ues of the fields <£>o an d $ a are obtained from the stationary-point conditions: 
= and = 0. The first condition is trivially fulfilled for $ a = 0, and 
the second one has a nontrivial solution $ = 5*0, which satisfies the equation 

a 2 = 2v [g (S ) + (5b)) , (5) 

with g (So) and g' (So) given by the following integrals: 

a' (S ) - f ^ R(k2) (7) 

Above, we have used the notation 

v = 2N c N f , D(k 2 ) = k 2 + M 2 (k 2 ), 

M(k 2 ) = R(k 2 )S + m, R(k 2 ) = r 2 (k 2 ), (8) 



where A^ c = 3 is the number of colors and Nf = 2 is the number of flavors. 

irk cone 
explicitly gives 

r d A k M(k 2 ) 



The quark condensate (for a single flavor) is defined as (qq) = ^-j^ , which 
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If m ^ 0, the above expression is ultraviolet divergent. However, in this case 
one should subtract the perturbative vacuum part from @, i.e. the piece with 
So = 0. The results is well-defined: 



{qq) > (qq) - (qq) s Q = -4JV, 



d A k ( M(k 2 ) m 



' I (2tt) 4 \k 2 + M 2 (k 2 ) k- ■ ///- 
d A k {2m + R{k 2 )S )R{k 2 )S 



-4N c S 09 '(S ) + 4N cmj (2?r)4 D{k2){k2 + m2) 
= -4N c S g'(So) + 0{m) (10) 

The model has three parameters: the cut-off A, the current quark mass m 
and the coupling constant -j. Two conditions for parameters are obtained by 
demanding that the pion decay constant, F n , and the pion mass, m^, assume 
their physical values: F n = 93MeV, and = 139MeV. To leading order in 
m the expressions for these quantities in our model are 



2 2uS mg'(S ) 

m - = p2 ' ( 12 ) 

We note that the elimination of g'{So) from r.h.s of the expression (|T^) , with 
help of the equation (|T0|), leads to 

2 m(qq) 

m - = ^2—' ( 13 ) 

in which we recognize as the Gell-Mann-Oakes-Renner relation. One parame- 
ter of the model is left free. Through the use of the stationary-point condition 
we can trade it for the value of So, which now is treated as the only remain- 
ing parameter. In principle we could still fit it, making use another helpful 
quantity, namely the gluon condensate. In our model it can be written as p( 



^^) = ^a 2 S 2 - m (qq). (14) 

However, the dependence of the gluon condensate on S (with F n and m w fixed 
at physical values) is very weak, which makes the fit impossible in practice 
due to a large uncertainty in our knowledge of (a s /ii G^G^). Reasonable 
parameters yield the value of the quark condensate (qq) ~ (250MeV) 3 , the 
current quark mass m ~ 5MeV and the gluon condensate (a s /n G^G^) ~ 
(330 ± 30MeV) 4 . This favors S /A in the range 0.2 - 0.5. 

Now we pass to the analysis of the temperature dependence of the quark 
condensate. For this purpose we use the Matsubara formalism, which leads to 
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the following simple thermalization rule for one-quark-loop integrals 



d A k [dk f d 3 k ~ f d 3 k 



where T is the temperature, / is the integrand of the quark loop, and k 
denotes the momentum running around the loop. The summation runs over 
the appropriate fermionic Matsubara frequencies, Ej = (2j + 1)ttT. After 
thermalization and integration of the angular dependence, the equations (|6]) 
and (|7] ) acquire the following finite-temperature form: 

»™ = 7 T ±h »™!M { e + w (16) 



Of course the value of So depends on T via the finite-temperature stationary- 
point condition: 



(Tft \ 
9t (Sq) + —g' T (So) J 



where now k denotes the length of the 3-momentum k. The finite-temperature 
expression for the quark condensate is 

(qq) T = -AN c Sog' T (So)+ (19) 
4iV c m - r k 2 {2m + R(k 2 + El)So)R(k 2 + El)So 



-T X; fdk 



it 2 ^ J (k 2 + El + M(k 2 + El) 2 )(k 2 + E 2 +m 2 



There is an important remark to be made in this place. The described way of 
thermalization, typically done in effective chiral quark models, tacitly assumes 
that the model parameters (here a and A) do not explicitly depend on T. 
There is no sound reason for this assumption in our effective theory, except for 
simplicity. Allowing temperature dependence of model parameters influences 
the values of T c . We will return to this question shortly. 

Table 1 shows our numerical results, with the meaning of various quantities 
given in the caption. The striking feature here is the very weak dependence 
of the critical temperature T c on the model parameters and on the regulator 
used. We find T c ~ lOOMeV, which is somewhat less than in conventional 
Nambu-Jona-Lasinio model, where it is of the order of 120 — 150MeV. The 
reason for this lower value is the following. The integrand of the expression 
for gx in local models is simply k 2 /(k 2 + E 2 + M 2 ), where M is the quark 
mass, instead of our expression k 2 R 2 (k 2 + E 2 n )j{k 2 + E\ + M 2 R 2 (k 2 + E 2 )). 
The presence of the regulator function R causes the much more rapid decrease 



5 



Table 1 

Results for the two considered regulators and various parameter sets. The first 
column is the ratio So/A, the next three columns display the values of the model 
parameters A, a, and m, adjusted in such a way that F n = 93MeV and m n = 
139MeV for the given So/A. The next two columns show the calculated values of So 
and {—qq)3 at T = 0. The last three columns give the results obtained in the strict 

chiral limit, with m set to 0, and A and a from columns (2-3). The quantities Sq m= o 

i 

and {—qq)m=o denote the values of the scalar field and the quark condensate for this 
case and for T = 0, and T cr it denotes the critical temperature of the second-order 
chiral phase transition. 



So 
A 


A 


a 


m 


S 


(-<?<?> ^ 


5*0,m=0 


i 

(-qq)L=o 


-^crit 




MeV 


MeV 


MeV 


MeV 


MeV 


MeV 


MeV 


MeV 


Instanton-model re 


gulator 












0.2 


1506 


219 


3.5 


301 


287 


256 


274 


102 


0.3 


1191 


165 


4.9 


357 


256 


310 


246 


101 


0.4 


1019 


134 


6.2 


408 


238 


358 


229 


100 


0.5 


907 


115 


7.3 


454 


225 


401 


218 


99 


Gaussian regulator 


0.2 


1155 


296 


3.9 


231 


290 


192 


278 


106 


0.3 


894 


218 


5.9 


268 


250 


230 


243 


105 


0.4 


754 


176 


7.7 


302 


227 


263 


221 


104 


0.5 


666 


149 


9.4 


333 


211 


295 


207 


103 



of the integrand with the temperature T, which enters through the Matsubara 
energies E n . The more rapid decrease of g? with T leads directly to a lower 
critical temperature T c as compared to the local models. One should bare in 
mind here that any explicit dependence of a on T, not considered here, would 
influence the values of T c . For instance, an increase of a by 20% (equivalent 
to decreasing the coupling constant of the four-quark interaction) leads to 
an increase of T c to values around 120MeV. We note that calculations in 
instanton-based models of the QCD vacuum [|l]|,[n],[32|,[33|] yield values of T c in 
the range 125-175MeV. 



To summarize this part, we note that at the one-quark- loop level the behavior 
of the non-local model is qualitatively not different from the local models. We 
find the chiral phase transition, which is second-order in the strict chiral limit, 
and the corresponding critical temperature is very weakly dependent on the 
details of the model. 

In the remaining part of this paper we consider meson loops effects. Such 
effects, although suppressed by 1/N C , are expected to be important due to 
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Fig. 1. The temperature dependence of So and {—qq) 3 for the instanton-liquid- model 

regulator with a = 134MeV, A = 1019MeV and m = 6.19MeV or m = 0. The 

meaning if the curves is as follows: So - solid line, {—qq)z - long-dash line, So )m =o _ 

i 

medium-dash line, {—qq)m=o ~ short-dash line. Note considerable difference between 
the strict chiral limit and the fmite-m case. 



the low mass of the pion. We follow the method described in Refs. fl27| , |28|l , 
where the meson loops have been introduced in effective chiral models in a 
consistent (i.e. symmetry-conserving) way. As mentioned in the introduction, 
conventional Nambu-Jona-Lasinio models require an extra regulator for the 
meson loop - making the quark loop finite via momentum cut-off, the proper- 
time or Pauli-Villars method, still leaves the momentum in the meson loop 
unconstrained. With non-local regulators this is no longer case. The interac- 
tions @ cut-off any momentum flowing through a quark vertex. As a result, 
calculation of a diagram with any number of loops is finite. 



In particular, meson loops are finite, and we now incorporate them in our 
study. First we analyze the problem at T = 0. The meson-loop contribution 
in the effective action can be obtained in a standard way by integrating out 
meson fluctuations around the stationary-point configuration P8 |. The new 
stationary-point condition for Sq assumes the form 



X 




= 



Fig. 2. The diagrams contributing to the stationary-point equation ( |20| ) with meson 
loops included. The dashed line corresponds to the meson propagators K n and K a 
defined in Eq. (|2"3|). The cross denotes the term a 2 So. 
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S [a 2 - 2vg(S ) - 2iy^-g'(S ^ 



+ ( V *W K °W + n bVMKM) = , (20) 

where Nb = Nj — 1 = 3 is the number of pions, K c and K n are meson 
propagators at the quark-loop level and the 3-point quark loops (see Fig. 2) 
are equal to: 

V 'V> = - 2V I J^ D^)D\k + P ) { - 2M{k) k ' ik + P) 

-M(k + p)(k 2 + M 2 {k))) (21) 
d 4 k R 2 (k)R(k + p) 
(2tt) 4 D 2 {k)D{k + p) 
+ M(k + p)(k 2 -M 2 {k))) (22) 



T ,/x n r d A k R 2 (k)R(k + p) / „ 



The inverse meson propagators can be written as: 
X- 1 (p) = ^ 2 / 7r (p) + £2^ ), 

^ 1 b) = ^ 1 (p)+4</.(p). (23) 
where we have introduced 



2_ r d 4 k / R(k)R(k+p)(k-(k + p)+M(k)M(k + p)) 
J*W- p2 J /( 2 7r) 4 l v D(k)D(k + p) 

#m\ (24) 



and 



f(\ 1 /• rf 4 fc i?(A;)i?(£; + p).M(£;).M(A; + p) 

Mpj ~ 5 2 y (2vr) 4 D(k)D(k+ P ) ■ [b) 

All terms included in the stationary-point condition are shown in Fig. 2. A 
similar expression can be obtained for the quark condensate. 

The results are somewhat surprising. We have found large cancellations in the 
meson-loop contribution, such that the net correction to the value of So is less 
than 1%, and to (qq) of the order of 3-4% on top of the quark-loop value. To be 
more precise, for the Gaussian regulator with the parameters A = 754MeV, 
a = 176MeV, and m = we find (qq) = (0.0505 + 0.0021 - 0.0041)GeV 3 , 
where we have decomposed in the quark-loop, a-loop, and the pion-loop con- 
tribution, respectively. We note the cancellation between the a- and pion-loop 
contributions. The phenomenon occurs for various regulators and parameter 
sets. It is also instructive to look at the integrand of the p-integration in Eq. 
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(|20|). After carrying the angular integration, we are left with a radial integra- 
tion over \p\. The integrand in the pion-loop case changes sign as we increase 
\p\ , and this results in an additional cancellation. At large values of \p\ the in- 
tegrand drops to zero. This behavior is quite different from the calculations in 
the local model, as in Ref. [p5| . There the analogous integrands in the meson- 
loop contributions increase with \p\, and the results depend strongly on the 
value of the momentum cut-off in the meson loop. 

Having found small meson-loop effects for T = we expect them to remain 
small at finite T. This is indeed the case. The expression ( |20| ) and the similar 
expression for the quark condensate can be thermalized by the prescription 



dig) . The only difference now is that in the mesonic loop we have to use 



bosonic Matsubara frequencies, i.e. in the meson-loop we replace 

with ujj = 2njT. The numerical calculations show that the effect of the me- 
son loops for finite temperatures is tiny, and the quark loop is dominant. The 
critical temperatures are practically unchanged. This is different than in lo- 
cal models, where the meson loops have a substantial effect on the finite-T 
behavior of the model ||29[| . 



Our main conclusions are as follows: the restoration of chiral symmetry in 
non-local chiral quark models proceeds similarly to the conventional Nambu- 
Jona-Lasinio models. In the strict chiral limit the transition is second order, 
however the critical temperatures are lower, around lOOMeV. For finite current 
quark masses we find a smooth cross-over to the restored phase. An interesting 
result is a very small contribution of meson loops to the quark condensate, 
which is in contrast to the calculations within the Nambu-Jona-Lasinio model. 
The effect occurs for a variety of regulators used. It is encouraging, since it 
indicates that we can rely on the simple-to-do one-quark-loop calculations. It 
remains to be seen whether in non-local models the meson loop effects are also 
negligible for meson correlators and other Green's functions. 
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